Collective oscillations of an interacting trapped Fermi gas 
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We calculate the effects of two-body interactions on the 
low frequency oscillations of a normal Fermi gas confined in 
a harmonic trap. The mean field contribution to the collec- 
tive frequencies is evaluated in the collisionless regime using a 
sum rule approach. We also discuss the transition between the 
collisionless and hydrodynamic regime with special emphasis 
to the spin dipole mode in which two atomic clouds occupy- 
ing different spin states oscillate in opposite phase. The spin 
dipole mode is predicted to be overdamped in the hydrody- 
namic regime. The relaxation time is calculated as a function 
of temperature and the effects of Fermi statistics are explicitly 
pointed out. 
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The investigation of collective excitations in trapped 
atomic gases has become an active research field, stim- 
ulated by the experimental realization of Bose-Einstein 
condensation 0] . Because of the high density of the con- 
densate, interaction effects in cold Bose gases are crucial 
H in order to explain the experimental results for the 
collective frequencies [|]. In Fermi gases the Pauli exclu- 
sion principle makes the density of the trapped gas more 
dilute, thereby reducing the effects of interactions. The 
equilibrium properties of trapped Fermi gases as well as 
of Fermi-Bose mixtures have been already the object of 
theoretical calculations based on mean field approaches 
. On the other hand experiments aiming to produce 
samples of Fermi gases in conditions of quantum degen- 
eracy are also becoming available ||. 

In this letter we provide a theoretical discussion of in- 
teraction effects on the collective oscillations of a trapped 
Fermi gas. The very high precision of frequency measure- 
ments, the large value of the scattering length exhibited 
by some atomic species and the possibility of pointing 
out effects of quantum statistics in the collisional term, 
make the study of collective oscillations in these systems 
a promising area of investigation. We will consider Fermi 
gases occupying two distinct spin states (hereafter called, 
for simplicity, spin up and spin down respectively). In 
fact only in this case can the effects of the interaction 
generated by s-wave scattering be explored. We will fur- 
ther limit the discussion to the case of excitations of low 
multipolarity in the normal phase, corresponding to the 
easiest realization of future experiments. In fact the tran- 
sition to the superfluid phase is predicted ||,|l0| to occur 
at very low temperatures. For a discussion of the oscil- 
lations in the superfluid phase see . 

At very low temperatures collisions are quenched by 



Fermi statistics and the system is in the collisionless 
regime, corresponding to the propagation of zero sound 
in traditional Fermi liquids fT2| ], A useful approach to 
the study of collective excitations in this regime is pro- 
vided by sum rules p^ j. In this approach the excitation 
energies are estimated through the ratio 
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between the k = 3 and k = 1 moments mj. = 
J S(E)E k dE of the dynamic structure factor 

S(E) = £|<n|f|0)| a *(£-7iw n ) (2) 



relative to a given operator F. For simplicity we have 
considered the T = case, but the formalism can be 
naturally extended to finite temperature. The moments 
mi and m.3 are easily evaluated in terms of commutators. 
In fact, using the completeness relation, one can write 



mi =^(0\[F^[H,F]]\0), 
m ,=\{Q\[[F\H],[H, [H,F]}}\0), 
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thereby avoiding the most difficult problem of determin- 
ing the full function S(E). In eq.(3-4) |0) is the ground 
state of the many body system. Let us consider a sym- 
metric configuration with — N± — N/2 and the same 
confining potential Vj xt — V} xt = V ext for the two spin 
species. We will evaluate the commutators using the 
mean field hamiltonian 



H 



4^ 2m 



l) 



(5) 



where g = 4irh 2 a/m is the interaction coupling con- 
stant fixed by the s-wave scattering length a. By choos- 
ing an isotropic harmonic potential V ex t( r ) = h mu! ho r2 
the commutators with the monopole (F — J2i r i)i 
quadrupole (F = J2i r i ~ %zf) and dipole (F = J2i z i) 
operators can be easily evaluated and one obtains the 
following results (T3 
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LU D = LU ho (8) 

where Ei nt = g J p{r) 2 dx/A is the mean field interaction 
energy, Eho — mu} 2 Jr 2 p(r)dr/2 is the potential oscilla- 
tor energy, E kin = £\ (0|p?/2m|0) is the kinetic energy 
relative to the ground state and p = Pt + Pi is the total 
density of the gas. In deriving the second equalities in 
(fl@) we nave used the virial theorem Q 

2-Efcm — 2Eho + 3E int = (9) 

which allows one to calculate the deviations of the col- 
lective frequencies from the ideal gas prediction 2u>ho 
in terms of the ratio Ei nt /Eh - Differently from the 
quadrupole and the monopole, the dipole frequency is 
not affected by interactions. In fact this mode corre- 
sponds to the oscillation of the center of mass of the 
gas driven by the external harmonic potential. It is 
worth noticing that results (§|8|) hold also in the case of 
trapped bosons |H] , provided one uses the corresponding 
expression for the interaction energy. In the case of Bosc- 
Einstein condensed gases the ground state kinetic energy 
is strongly quenched by interactions and, for large TV, the 
ratio Eint/ Eho is equal to 2/3 (see eq. ([)))). In this limit 
one recovers the results ujm — V^^ho and ujq = \/2cjh 
predicted by the hydrodynamic theory of superfluids fll4[ . 
Conversely, in the Fermi case the interaction energy is, in 
most cases, only a small perturbation that can be safely 
estimated using the Thomas-Fermi expression 

p{r) = (1?) 3/2 i (Mo ~ y -*w) 3/2 ( 10 ) 

for the ground state density. In this equation [Iq = 
(SNy^huJho is the chemical potential of the Fermi gas, 
fixed by the normalization condition. Using expression 
( |l0| ) for the density one obtains the result 

E int NWa 

— - =a (11) 

J^ho aho 

where a = 8192 s/2 3 1 / 6 /2835tt 2 ~ 0.50 while a ho = 
(H/muJho) 1 ^ 2 is the usual harmonic oscillator lenght. 
Interaction effects are governed by the combination 
N x ^aj 'ah showing that in order to emphasize the role 
of interactions it is much more efficient to increase the 
ratio a/ aho rather than the number of atoms. The ratio 
(|ll|) can be also casted in the form Ei nt / E) lo — 0.3 x kpa 
where kp = {2ra{?>N) 1 /' i ijJho/K) 1 ^ 2 is the Fermi momen- 
tum of the trapped gas |l5| . 

In analogous way one can calculate the frequency of the 
out of phase oscillations, hereafter called spin excitations. 
We report here the result for the most relevant spin dipole 
mode excited by the operator F = Zi — z i- This 
mode corresponds to a relative oscillation of the centers 
of mass of the spin-up and spin-down clouds and is the 
analog of the giant dipole resonance exhibited by atomic 
nuclei p6|. The frequency of the spin dipole mode is 



found to be 

usd = ^"ho-^Jj J\d z p\ 2 dv 

~u ho \l-a' (12) 

V a ho J 

where a' = 128 %/2 3 1/6 /35tt 2 ~ 0.63 and, in the sec- 
ond equality, we have used the ground state density ( |l0| ) 
to evaluate the interaction contribution to first order in 
the scattering length. Notice that interactions affect the 
spin-dipole mode through the same combination of pa- 
rameters characterizing the ratio ( [ill) • Typical values for 
the fermionic isotope of potassium show that interaction 
effects are rather small. For example, using N — 5 x 10 5 
and a I aho = 6x 10~ 3 one finds that the dipole frequency 
decreases by ~ 3%. Significantly larger corrections are 
predicted in the case of 6 Li where the scattering length is 
a factor 10 larger and negative. Notice that the integrals 
characterizing the interaction contribution to the collec- 
tive frequencies, are calculated here at T = and would 
decrease at higher temperature. We also note that, at 
T = 0, sum rules provide only an upper bound to the 
frequency of the lowest state excited by the operator F. 
However, if one considers first order corrections in the in- 
teraction constant a, the sum rule bounds (^^,12) can be 
shown to correspond to the exact value of the collective 
frequencies. 

The results derived above can be easily generalized to 
the case of deformed traps in which the external potential 
reads V ext (r) — mui\ (x 2 + y 2 + X 2 z 2 )/2 and A = ll> z /uj± 
is the deformation parameter of the trap. By assum- 
ing that interaction effects are smaller than the unper- 
turbed splitting between the radial and axial frequen- 
cies, one finds simple results also in this case. The 
new decoupled frequencies, associated with the radial 
and axial excitation operators F ra di a i — X)j x f + Vi an d 
Faxiai = J2i z i> become, to first order in a, uj radla i = 2uj±_ 
and LUaxiai — 2oj z {\ — 3/16 aiV 1 / 6 a/ 'aho) showing that 
only the axial mode is affected by the interaction. The 
oscillator length a/ lo = (h/muj) 1 / 2 is here defined in 
terms of the geometrical average ZU = w^A 1 / 3 of the 
three frequencies. For the spin dipole mode result (12) 
is easily generalized to both the radial and axial direc- 
tions. One finds u> ra dial = — a'N^-^a/aho) and 

Uaxial = U>a(l - a! ' N 1 / 6 'a / ' a ho ) ■ 

So far we have ignored the effects of collisions. If the 
collisional frequency is much larger than the frequency 
of the collective excitations then the system is in the 
hydrodynamic regime, also known as first sound regime. 
An important question is whether the transition between 
the collisionless and hydrodynamic regimes takes place in 
the degenerate or classical regime for our trapped gases. 
Let us first discuss the collective frequencies of a trapped 
Fermi gas in the full hydrodynamic regime (see |l7| for a 
recent discussion). In the spherical case the quadrupole 
frequency becomes luq D = V2u>ho instead of (0). This 
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result is independent of statistics and holds also for a 
classical gas |Lq] . The reduction with respect to the col- 
lisionless value is due to the fact that, in the collisional 
regime, the only restoring force for surface excitations 
arises from the external field. For the monopole fre- 
quency one instead finds that result @ holds also in the 
hydrodynamic regime. The above analysis can be also 
extended to the case of deformed traps where one 
finds that, in the absence of mean field effects, the col- 
lective frequencies coincide with the ones holding for a 
classical gas Q . 

Collisions are expected to have more dramatic conse- 
quences on the spin-dipole oscillation since they do not 
conserve the spin current and consequently give rise, in 
the hydrodynamic regime, to a pure diffusive mode. The 
investigation of the spin dipole mode is consequently ex- 
pected to be a sensitive test of the role of collisions and, 
possibly, of quantum statistics. The equations for the 
spin dipole oscillation can be easily obtained starting 
from the Boltzmann equation where the effects of Fermi 
statistics are included in the collisional integral. For a 
first estimate let us ignore the mean field effect, which 
is responsible for the frequency shift of (12). Using the 
method of the averages recently developed in |2(J for clas- 
sical trapped gases, one obtains the following coupled 
equations: 

3t(*T _ z i) ~ ( v *1 ~ v zl) = ( 13 ) 
d t (v z f - v zi ) + u 2 z { Z] - Z[ ) = {{v z] - v zi )I coU ) 

(14) 

where the average (...) is taken both in coordinate and 
momentum space. The collisional term is given by 

((i>zT - v zi )I col i) = 
am 6 2 f 

- J^eJJ / drdv 1 dv 2 dtl\v 1 - v 2 \(v lz - v 2z ) x 

x [(1 - / 1T )(1 - / 2i )/( T /^ - / 1T / 2i (l - / 1T )(1 - f 2l )\ 

(15) 

where a = Aira 2 is the total cross section, Vj and v' 2 are 
the velocities of the particles 1 and 2 after the collision, 
and /f, fi are the distribution functions relative to the 
two spin components, normalized to (m/h) 3 J f^dvdv — 
(m/h) 3 J fidrdv = N/2. The collisional term can be 
estimated by assuming that during the oscillation the 
distribution functions of the two spin species behaves, 
in velocity space, as f n (v Xl v y ,v z ) = fo(v x ,v y ,v z ± u), 
where /o is the distribution function of each compo- 
nent at thermal equilibrium. This corresponds to a rigid 
displacement of the velocity Fermi distributions of the 
two spin components in opposite directions. One finds 
(w z | — v z i) — 2u and, by developing the integral (15) to 
first order in it, one can finally write 



((VzT ~ Vzl)Icoll) = - 
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where 
1 



drdvidv 2 dn\v\(v — v 
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t 24irh e NK B T , 

x /o(vi)/o(v 2 )(l - /„K))(1 - /oK)) (17) 

defines the relevant relaxation time of the spin dipole 
oscillation. In (17) we have defined v = vi — v 2 and 
v' = Vj — v 2 . As a consequence of ( fl6|) the equations 
of motion (13-14) for the spin dipole oscillation take the 
simple form of a damped harmonic oscillator. Looking 
for solutions of the form e~ %U3t the dispersion law is given 
by 
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(18) 



showing that the oscillations become overdamped if 

UJ Z T < 1/2. 

The relaxation time (17) is easily evaluated at high 
temperature where the effects of Fermi statistics are neg- 
ligible. In this case it takes the form 
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(19) 



where j c i — vthcrp(0)/2 is the classical collisional rate, 
v th = (SKbT/ttiti) 1 ^ 2 is the thermal velocity and p(0) = 
(0) + pi(0) is the central density of the gas. At tem- 
peratures smaller than the Fermi temperature Tp — 
(SNy^HuJho/ Kb the effects of statistics become impor- 
tant and, for T — > 0, the relaxation time becomes larger 
and larger exhibiting the typical behaviour of Fermi sys- 
tems. A similar behaviour has been recently shown to oc- 
cur in the relaxation of the motion of a classical particle 
inside a degenerate trapped Fermi gas (l9) . By a proper 
change of the variables entering the collisional integral 
(17), the dimensionless quantity uj z t can be written in 
the useful form (for oscillations along the radial direc- 
tion one should consider the quantity ujj_t = uj z t/X) 
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UJ Z T 



_±_ A -2/3 ^1/3^ 

3 4 / 3 7r 1 a ho 



F 



T 



(20) 



where F(t) is a dimensionless function, sensitive to Fermi 
statistics, which determines the temperature dependence 
of the relaxation time. This function is plotted in Fig.(Q). 
For large values of the reduced temperature t = T/Tp it 
approaches the classical behaviour 1/t, while for t smaller 
than 1 one observes, as expected, important deviations 
due to quantum effects. At very low temperatures the 
function F(t) can also be calculated analytically and we 
find 



2.2 



F{t) = 8ir 2 t 



(21) 



Of course result (|21j) holds above the superfluid transi- 
tion. The figure shows that the system will be in the col- 
lisionless regime both at sufficiently high and low temper- 
atures. In the first case the gas is classical and collisions 
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are rare because the density is very low. In the second 
one the gas is degenerate and collisions are rare because 
of Fermi statistics. The possibility for the system to reach 
the overdamped regime lu z t < 1/2 depends in a crucial 
way on the value of the physical parameters entering (20). 
The maximum value of F(t), reached at t ~ 0.4, is ~ 1.6 
so that the condition for the existence of overdamping 
is given by A 1/ ' 3 7V 1 / 3 |a|/a; lo > 2.1 (for the spin dipole 
oscillation along the radial axis the condition is instead 
\-i/e N i/3\ a y aho > 2.1). This condition should be easily 
achievable in the case of lithium. For example, choosing 
a spherical trap with N = 5x 10 5 and a/a^o = —3 x 10 -2 
one finds overdamping for 0.2 < T/Tp < 0.8. For lower 
temperature the system will exhibit a transition to the 
collisionless regime where the spin oscillation is damped 
according to eq. (|l8|) . In the case of 40 K, where the value 
of the scattering length is much smaller, the system is ex- 
pected to be always far from the hydrodynamic regime. 
Nevertheless an accurate investigation of the damping 
might point out the effects of Fermi statistics also in this 
case. In fact the damping of the spin oscillation is pre- 
dicted to strongly decrease at low temperature in con- 
trast to the classical behaviour. Notice that the above 
estimates are based on the unperturbed value lj z (u>±) 
for the spin dipole frequency. At low temperature one 
can easily renormalize these values by including the mean 
field effect through eq.(12). 

In conclusion, we have provided a systematic discus- 
sion of interactions effects on the propagation of collective 
oscillations in trapped Fermi gases. In the collisionless 
regime interactions renormalize the collective frequencies 
through mean field effects. The predicted shifts should 
be visible in experiments due to the high precision of fre- 
quency measurements. Collisions are responsible for the 
damping of the oscillations and their effect turns out to 
be very sensitive to Fermi statistics. Overdamped oscilla- 
tions are predicted to occur in the case of the spin dipole 
oscillation and should be visible in lithium due to the 
high value of the scattering lenght. The present analysis 
can be naturally extended to include asymmetric config- 
urations with different numbers of atoms in the two spin 
states (iVf 7^ iVjJ and/or different trapping potentials 

(vL * vit). 

Useful discussions with D. Guery-Odelin, H. Stoof and 
F. Zambelli are acknowledged. 
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FIG. 1. Spin dipole relaxation function (see eq.(|20[)) as a 
function of the reduced temperature t = T/Tf- The classical 
prediction 1/t (dahsed line) is also shown. 
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